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. Abstract 
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The structure of the general solution of the inhomogeneous evolution equa- 
tions allows the separation of a photon structure function into perturbative 
(JL, ■ ( "anomalous" ) and non-perturbative contributions. The former part is fully 

calculable, and can be identified with the high-mass contributions to the 
dispersion integral in the photon mass. Properly normalized "state" distri- 
butions can be defined, where the 7 — > qq splitting probability is factored 
out. These state distributions are shown to be useful in the description of 
the hadronic event properties, and necessary for a proper eikonalization of 
jet cross sections. Convenient parametrizations are provided both for the 
state and for the full anomalous parton distributions. The non-perturbative 
parts of the parton distribution functions of the photon are identified with 
the low-mass contributions to the dispersion integral. Their normalizations, 
as well as the value of the scale Qq at which the perturbative parts vanish, 
are fixed by approximating the low-mass contributions by a discrete, finite 
sum of vector mesons. The shapes of these hadronic distributions are fit- 
ted to the available data on F%{x, Q 2 ). Parametrizations are provided for 
Qq = 0.6 GeV and Qq = 2 GeV, both in the DIS and the MS factorization 
schemes. The full parametrizations are extended towards virtual photons. 
Finally, the often-used "FKP-plus-TPC/27" solution for F^{x, Q 2 ) is com- 
mented upon. 



Heisenberg Fellow. 



CERN-TH/95-62 



1 Introduction 



Inclusive cross sections, such as deep-inelastic electron-photon scattering e^y — > e + X 
or high-p^ jet production in two-photon collisions 77 — > jet + X, are fully expressed in 
terms of the parton distribution functions (PDFs) of the photon f2{x,Q 2 ) (for recent 
reviews on the photon structure functions see e.g. refs. Jj], ^, p[ f|, fjj, ||). Correspondingly, 
such measurements allow, in principle, a unique determination of the PDFs of the photon. 
However, already for the simplest determination, namely the extraction of F^x, Q 2 ) from 
deep-inelastic electron-photon scattering, additional information about the hadronic final 
state is necessary ||. The reason is simple: in contrast to lepton-nucleon scattering, 
where (x, Q 2 ) are fixed by the measurement of the outgoing lepton solely, in ej scattering 
x usually can only be determined through a reconstruction of the hadronic energy W 
from the visible one by virtue of the relation x = Q 2 /(W 2 + Q 2 ). Clearly, an exclusive 
description of the event properties is not only of interest for acceptance corrections but 
provides also additional insight in the structure of the photon. 

In this paper we shall investigate the photon-to-qq splitting probability and define 
(normalized) "state" distributions where this probability factor is split off. We shall 
show that the use of these state distributions is of great advantage in obtaining a correct 
description of the hadronic final state, notably of the photon remnants and the initial- 
state gluon radiation in parton-shower programs. We shall see that the state distributions 
enter also the calculation of jet cross sections at high energies where unitarity corrections 
have to be taken into account. 

Furthermore, we shall show how to obtain a constrained parametrization of the (in- 
clusive) PDFs of the photon. Current parametrizations either are fully unconstrained, i.e. 
all parameters are fitted to data 0, or rely on severe prejudices about the input dis- 
tributions, such as approximating these by the ones of the pion || [H| or modelling them 



by the quark-parton results (box-diagram) with the quark masses as free parameters JTT 
By contrast, the shapes of our input distributions are determined by the K7 data, while 
the normalizations are constrained. There exists a one-to-one relation between the size of 
the non-perturbative part and the scale Qo introduced to separate the non-perturbative 
and the perturbative parts. This scale should be universal, i.e. process-independent and 
can, for example, be constrained by a study of the 7p total cross section We also 



investigate the sensitivity of the PDFs to the photonic factorization scheme associated 
with the presence of a direct contribution C 7 to F^- The current analysis updates our 



old parametrization of the photon PDFs [12 



Next we discuss the dependence of the PDFs of the photon on its mass (virtuality) 
P 2 , i.e. the so-called target-mass effects. Our parametrization can, in fact, also be used 
for the PDFs of a virtual photon, and we shall outline the assumptions that are behind 
this. Effects of a non-zero virtuality have attracted considerable interest [|l^, [14], f|, |15|, 



TBI, O, ITS], but no explicit (analytic) parametrization was available up to now. 



The outline of the paper is as follows. In the next section we first discuss the photon- 
to-qq splitting and the issues of introducing a scale Q to separate low-mass contributions 
from perturbatively calculable high-mass contributions. The latter, the so-called anoma- 
lous contributions, are then expressed in terms of "state distributions" that enter the 
description of the hadronic event (section 2.1), the eikonalization of jet cross sections 
(section 2.2), and the virtual PDFs (section 2.3). Section 3 contains details about the 
parametrizations of the parameter-free state distributions and the anomalous distribu- 
tions. The small- 2 behaviour of the distribution functions is outlined. In section 4 we 
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discuss the normalization of the non-perturbative (low-mass) contributions to the pho- 
tonic PDFs ("hadronic distributions") and compare our fits to the data. Finally we 
comment upon the often-used approximation of K7 by the sum of the (Q 2 -independent) 
TPC/27 parametrization [19) and the QCD expression of FKP |20 |. Our results are sum- 



marized in section 5. There we also discuss further aspects of our separation of the photon 
distribution functions into long-distance and short- distance components, in particular the 
connection with the representation of a photonic PDF as a dispersion integral in the 
photon mass. 



2 Separation of the non-perturbative part 

Perturbative QCD predicts only the Q 2 evolution of the PDFs of the photon fi(x,Q 2 ), 
via a set of inhomogeneous differential equations of the first kind. Hence the solutions 
fi(x, Q 2 ) require the specification of the PDFs at some Q 2 = Qq. It is well known 
that these "input" distributions fi(x, Qq) have to have a considerable "hard" component 
if Qq is chosen to be a typical hadronic input scale (Q ~ 2GeV). This is quite in 
contrast to hadronic PDFs which vanish in a powerlike way oc (1 - x) p (p ~ 1) for 
x — > 1. Correspondingly, there is no unique guiding principle as to what the photonic 
input distributions should look like. Rather, parametrizations starting at large Q H [LlJ 
are obtained by just fitting the parameters of the input distributions f2(x, Qq) (shapes and 
normalizations) to the available F2 data. Moreover, since these data are so far restricted 
to large x, basically no constraint on the gluon distribution of the photon exists today in 
such an approach. 

Lowering the starting scale Qo to about 0.5-0.7 GeV, it is, however, possible to describe 
the F2 data with photonic PDFs obtained from a purely hadronic input []9], [12], [TO]. The 
latter is usually estimated by the PDFs of the pion using vector-meson-dominance (VMD) 
and the additive quark model. The recent next-to-leading-order parametrizations AFG 
pl and GRV M essentially differ in the ansatz to take the photon to be a coherent |M or 



incoherent M superposition of vector mesons. Both are parametrizations of the inclusive 



PDFs of the photon. Our aim [0 is to obtain PDFs where the probability for the photon 
to fluctuate into a hadronic state is made explicit. Firstly, this allows us to constrain the 
normalizations of the non-perturbative input distributions, while the shapes are fitted to 
F2 data. Secondly, we can trace the complete evolution of the parton showering, otherwise 
hidden in the expression of fi(x, Q 2 ). Thirdly, differences in the remnant structure of the 
photon can be respected. Finally, our parametrizations hold also for a virtual photon, to 
which they can easily be applied because the parametrizations are analytic in both P 2 
and Q 2 . On the other hand, our parametrizations are leading-order ones only. 
Analogously to our decomposition of the 7p total cross section 

7P _ 7P , 7P , 7P fl ^ 

°tot — °dir ' °VMD ' °anom > \ L I 

we decompose the PDFs of the photon as 

P a (x, Q 2 ) = f2^(x, Q 2 ) + p a ' YMD (x, Q 2 , Ql) + /r nom (z, Q 2 ; Ql) ■ (2) 

The ansatz ([[]), or eq. (@), is based on the assumption that there exists a scale Q which 
divides the spectrum of 7 — > qq fluctuations: above Qo these can be described pertur- 
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bativelyQ ("anomalous" term in eqs. ([!]) and (0)), while below Q the fluctuations are 
assumed to give vector-meson states. 

A reminder of the familiar scale ambiguity problem may be in place. In this paper 
we work throughout with the deep-inelastic-scattering convention in mind, where the 
Q 2 scale is set by the virtuality of the probing photon. Fluctuations 7 «-> qq of the 
probed photon are resolved as long as the "scale" k of these fluctuations is below Q. 
However, the precise physics definition of k is left unspecified. When PDFs extracted 
from deep inelastic scattering are applied to other processes, such as qq' — > qq' high-p^ 
jet production, the choice of relevant scale Q is ambiguous. It could be associated with the 
transverse momentum p± of the jets, or some multiple thereof. Conversely, our transverse 
momentum cut-off, used in |12| to separate between direct, simple VMD (i.e. p°, u, and 
4>) and anomalous processes, is of the order of the Q ~ 0.6 GeV used in this paper, but 
the two could well differ by some amount. 

The first term in eq. (|2p describes the (properly normalized) probability distribution 
of a photon to remain a photon 

r a d "(x,Q 2 ) = Z 3 5 a ,5(l-x) , (3) 

where 

Z 3 = l~ E / x {P a ' VMD (x, Q\ Ql) + P a '—(x, Q 2 ; Ql)} • (4) 

Properly speaking, also fluctuations 7 <-> £ = e, \i or r, should be taken into 

account for Z%. This contribution is fully perturbatively calculable. In practice, Z3 « 1 
is a sufficiently good approximation for all applications. 

The PDFs of the photon, being the solution of an inhomogeneous evolution equation 
(or, more precisely, a system of equations) can always be written as the sum of two terms 

Hi*, Q 2 ) ~ ft^ix, Q 2 ) = Pa' PT (^ Q 2 , Ql) + Pa' NP & ^ Ql) • (5) 

The first term is a particular solution of the inhomogeneous equation with the boundary 
condition 

P a PT (x } Ql,Ql) = 0. (6) 

The second term, a general solution of the corresponding homogeneous evolution equation 
needs a (non-perturbative) input distribution at Q 2 = Ql 

f^ P (x,QlQl) = F? P (x). (7) 

While mathematically the decomposition (|5|) is valid for any Ql, for physics reasons 
we want to identify fa' PT with the perturbatively calculable distributions f2 ,anom arising 
from the point-like coupling of the photon to a quark-ant iquark pair. These distributions 
can be expressed as an integral of "state" distribution functions fa' q ^{.x,Q 2 ; k 2 ): 

P—(x, Q 2 ; Ql) = ^ E K |r Q 2 ' k ") > ( 8 ) 

which obey the standard, homogeneous evolution equations with the boundary condition 
f2**(x, k 2 ; k 2 ) = f2^(x) = | (x 2 + (1 - x) 2 ) (5 aq + 5 m ) . (9) 



1 A similar definition of the perturbative part of F^x^Q 2 ) is due to FKP (£0|. A comparison is 
postponed to section 4.5. 
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The solution (|8|) can be understood as follows. The probability for the photon to 
branch into a qq state at some high scale k is given by (a em /27r) 2e^dk 2 /k 2 . Once the 
photon has fluctuated, the possibility of additional QCD evolution between k and Q is just 
given by the standard evolution equations. This evolution conserves total momentum, so 
that 

E f'dxxU^Q 2 ; k 2 ) = 1 (for Q > k) , (10) 

where the sum now runs over the gluon and all the quark and antiquark species. 

Similarly, f2' NP should describe the partonic content of the photon that has fluctuated 
to a vector meson, i.e. 

Pa YMD (x, Q\ Ql) = E ^Pa V ^ Q\ Ql) ■ (11) 

v Jv 

Here 47ra em / f v gives the probability for the photon to fluctuate into the vector meson V, 
while f2' v obeys a momentum sum rule just as f'2' qq in eq. (|T0|) . Normally we let the sum 
run over p°, uo and 0, with the first responsible for the bulk of the contribution. 
In passing we note that eq. (|i~0"D yields 

^3 = l-E^-^E2e?lnfg'| . (12) 



v 



fy 27T ^ q WO 



It is convenient to think of the anomalous component as a continuous spectrum of states, 
characterized by the flavour q and the scale of creation k, just as the VMD components 
give a discrete spectrum of states. The probability (a em /27r) 2e 2 dk 2 /k 2 in eq. (ffl) is 



therefore the equivalent of the VMD couplings 47ra em / fy in eq. (|TT|). 
2.1 Hadronic event properties 

For inclusive quantities, such as the high-p^ jet rate in 77 collisions, only the full parton 
distributions f2{x,Q 2 ) are of interest. However, for an exclusive description of event 
properties, states with different initial flavours and k values behave differently |21], [12]. A 



physically transparent description of the complete hadronic final state therefore requires 
that f2(x,Q 2 ) be available subdivided into the f2'™(x,Q 2 ) and / a 7 ' VMD (x, Q 2 ) parts, 
with these in their turn subdivided as an integral or sum of state distributions f2' q ^(x, Q 2 ) 
and f2' v (x,Q 2 ). Let us elaborate on this. 

1. Hard processes (at a scale p± ~ Q) in photon-induced reactions are given by the 
usual 2^2 partonic processes, where the parton a in the photon is selected according 
to the inclusive distributions f2(x,Q 2 ). In addition there are the processes where the 
photon interacts directly, i.e. via /J' dir (x, Q 2 ). For consistency, the direct cross sections 
have to be cut off at the same Qo as is used to define the anomalous PDFs, i.e. the cut-off 
will, in general, be different for different parametrizations and should be known! 

2. The 2 — > 2 partonic processes are then supplemented by initial (and final) state 
parton showers. In the case of partons in the photon, the inhomogeneous evolution equa- 
tions have to be used. However, parametrizations of the inclusive PDFs of the photon 
are not guaranteed to be positive-definite at the infrared cutoff used in the parton-shower 
programmes. Nor is the shape at Qo constrained in any way to correspond to a hadronic 
distribution, so that it makes sense to approximate the photon remnants by those of a 
hadron. On the other hand, both demands are fulfilled if f2 (x, Q 2 ) is split as in eqs. (H), 
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(§) and (|TT|). Moreover, existing parton shower programmes for homogeneous evolution, 
made and tested for hadr on- induced reactions, can then be used, since both f2' qq ~(x, Q 2 ; k 2 ) 
and f2' V (x,Q 2 ) evolve according to the homogeneous equations. The separation (0) also 
ensures that the "intrinsic" transverse momentum kj_ of the initial parton in the pho- 
ton can be correctly generated: a roughly Gaussian primordial k± distribution of width 
~ 0.5 GeV for low-mass hadronic fluctuations of the photon, and a dk\jk\ distribution 
for high-mass fluctuations. 



2.2 Eikonalization of jet cross-sections 



The state distributions f^ ,qq (x,Q 2 ; k 2 ) and f^' v (x,Q 2 ) enter also the calculation of total 
photon-induced cross sections at high energies. There (mini)jet cross sections rise much 
faster than allowed by the Froissart bound and unitarization corrections have to be taken 
into account. Unitarization effects can be estimated by the eikonal formula. Respecting 
the probability for the photon to be in different hadronic fluctuations, the correct formula 
for e.g. 7p interactions is 



a 



7P / 
inel \ 



V 

a 



Pv 



•cm 



2tt 



E2e 



dk 



d 2 b jl — exp 

d 2 b jl — exp 



2lm X Vp (s,b)}} 



-2Im X qm> (s,b,kl)]} ■ (13) 



The state distributions enter the hard cross sections, e.g. 



n- qqp (V k, 



do",- 



(!,, / dx 2 I dpif^(x 1 ,pl;k ± )ff(x 2 ,pi)^9(p ± -k 



(14) 



which in turn enter the eikonals x % m e Q- (0), e.g. x qqp — A mp a mp , where A qqp (b, k 2 ^) is 
an overlap function and a qqp (s, k\) = (x s q J(s, k\) + a qq r p d (s, k\). 



2.3 PDFs of a virtual photon 

Finally, the separation (0) is useful to study the effects of a non-zero photon virtuality 
P 2 . Generalized vector- meson dominance suggests the following dispersion relation in k 2 : 

The integration from zero to Q\ can be approximated by the (low-mass) vector-meson 
contributions, which show the usual fast P 2 fall-off predicted by pole dominance: 

/r' VMD (*,Q 2 ,P 2 ) = £ (-^) 2 ^=L f2> v (x,Q 2 ,P$) . (16) 

V \ m V ' / Jv 

The pole-dominance factor takes the reduced probability for a virtual photon to fluctuate 
into a vector meson into account. On the other hand, the PDFs of a virtual vector meson 
are completely unknown. The simplest choice, namely to identify them with those of 
a real pion, i.e. taking Pq = Ql in eq. ([16]) fll8| , leads to unphysically large hadronic 
contributions at small x. Moreover, this choice does not guarantee the vanishing of the 
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PDFs as P 2 — > Q 2 , as should be the case up to power-suppressed terms. Our choice 
Pq = max(P 2 , Qq) is based both on this limit and on the expectation that the PDFs 
should become more valence-like as P 2 increases. 

The anomalous contribution is given by the integration k 2 > Qq in eq. ( |i"5|) and we 
approximate it by 

(x, Q 2 , P 2 ) = f ^ ^ £ 2e 2 f2<*{x, Q 2 ; k 2 ) 

= f2>"(x,Q 2 } P 2 ) . (17) 
The ansatz ( |TTD guarantees that the virtual PDFs 

fX(x, Q 2 , P 2 ) = P a ' dil (x, Q 2 , P 2 ) + fT> YMD (x, Q 2 , P 2 ) + fT>"(x, Q 2 - P 2 ) (18) 
are correct in the region A 2 C P 2 C Q 2 : there the virtual PDFs can be calculated exactly 



and QCD predicts them to be given by eq. (ISI) with Qq replaced by P 2 |14{|. The PDFs 



of eq. (|T8| ) are exact also in the limit P 2 —>■ where they approach the real PDFs. 

On the other hand, for P 2 of the order of Qq, i.e. in the interesting "cross-over" region 
P 2 ~ m 2 , the P 2 dependence of the anomalous contribution is less certain. A power-like 
dependence might be present that could be estimated via eq. (|l|). Our choice of Pq 
in eq. (|T7D is motivated by (i) simplicity and (ii) the demand for a smooth transition 
towards higher P 2 , but also alternatives like P 2 = P 2 + Qq would have fulfilled these 
demands. Uncertainties arising from this choice may be estimated by varying P 2 , say, 
from max(P 2 , Q 2 /2) to max(P 2 , 2Q 2 ). 

As P 2 approaches Q 2 , the concept of virtual PDFs breaks down since powerlike cor- 
rections cx (P 2 /Q 2 ) p then become important besides the logarithmic ones. Nevertheless, 
the logarithmically enhanced terms of the P 2 — > Q 2 limit exhibit the correct behaviour: 
the quark distributions approach the box-diagram expression (quark-parton model re- 
sult) fL3| , and the gluon distribution vanishes faster than the quark ones Q. We finally 
want to emphasize that our parametrization of the virtual PDFs is analytic in all vari- 
ables and can hence be used easily. It is well known that hadronic distribution functions, 
such as f^ ,v (x,Q 2 , Pq), depend on the two momentum scales Q 2 and Pq only through 
the logarithmic integration of the coupling constant (see eq. (|19D below). However, also 
the solution of the inhomogeneous equation can be written in a similar form, see eq. 
below. 



3 The perturbative part 
3.1 The state distributions 

The distributions f2 ,q ^(x, Q 2 \ k 2 ) depend on two momentum scales, Q and k. However, 
the amount of evolution that occurs between these two scales is entirely characterized by 
the logarithmic integration of the strong coupling constant 

In leading order s = (1/6) ln[ln(Q 2 /A 2 )/ ln(fc 2 /A 2 )], where b = (33 - 2n f )/6. We have 
here assumed that the number of flavours is fixed, e.g. rif = 4, but it is straightforward to 
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split the Q 2 range into subranges with different rif values and correspondingly matched 
A.( n f) values. The number of flavours also enters in the fraction of the momentum carried 
by gluons to that carried by quarks; here a matching between different rif values is less 
transparent, but can still be attempted. 

The parametrization of the quark distributions is divided into a valence and a sea 
part. The sea part corresponds to the distributions /J' qq for o ^ q,q, while the valence 
part is the additional contribution obtained for a = q or a = q, i.e. 

fr,<n — f 7>qq _ /7>qq i ^7>qq 

Jq Jq Jq.val ' Jq,sea 

/q 7,qq = = f^m for q '_, q . (20) 

As always in QCD, the solutions are too complex to be given in closed form. One 
therefore has to resort to approximate parametrizations. We perform an evolution in x 
space choosing a strategy similar to the one proposed by Odorico P^| , wherein a large 
number of evolution histories are traced by Monte Carlo methods. The evolved parton 
densities are binned in x for several Q 2 values. Thereafter the x shape is parametrized in 
some simple form, and the Q 2 dependence of these parameters is in turn parametrized. 
The Minuit program [^3j is used to find suitable parametrization coefficients, but with 



frequent manual interaction to ensure a sensible behaviour. The fits have been made in 
the range 10" 4 < x < 1 and Q < Q < 2000Q (where Q « 0.5 GeV), but the forms 
have been chosen so that they can be used for all x. In regions where the distributions 
are large, the typical accuracy of the parametrizations is 1 %-2%. 
The parametrizations are 



_val „ ^val 



X&(X, S) = (cf X 2 + Cf (1 - Xf + Cf x{\ - X)) X C 4 (1 - X 2 ) C 5 

SGcl SGcl — , 

xf:; s f a (x,s)=crx c i (l-x) c s xf2™ l0 (x) 

xf^(x,s)=cfx c " ((l-x)(-\nx)) c hf Kl f(x) , (21) 



where 



xj:,:Zj-r) = 8 73x + 62x 2 (i - x )+(^-3 ] j x\nx + (2x - l)a;ln 2 x (22) 

_ A _|_ 7 r _|_ A r 2 

xf7£ q {x) = (1 - x) + 2x(l + x) lnx (23) 



and the coefficients q are given in table 3.1 



Note that the shapes and normalizations of the parametrizations become exact in the 
limit s — > 0: the expression in eq. @ is recovered in the limit s — > of the valence 
distribution, and the gluon (sea) distribution approaches eq. (|23|) (eq. (|22D ) times 4s (s 2 ). 
Hence, the sea distribution vanishes faster than the gluon distribution, which in turn 
vanishes faster than the valence one Q. This is because 4sfg'^ q (x) is obtained if the initial 
quark-distribution in eq. (|9]) is convoluted with the q — >• qg splitting kernel. Similarly, 
s Va sea,io( x ) i s obtained if the first-order gluon-distribution in eq. (|23|) is convoluted with 
the g — > qq splitting kernel. Moreover, our parametrizations exhibit the correct large-x 
behaviour, (1 — x) 8S//3 for the valence distribution, (1 — x) 1+8s ^ 3 for the gluon one, and 
(1 — x) 2+8s / 3 for the sea one. 

Analytic results can also be derived for the limit x — > 0. We find 

2s 



xfq^i(x, s) -f - exp ( — ) I (z v ) 



2 1 V 3 
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a 


b 


c 


d 


e 


/ 


9 


c l 


1.5 


— 


— 


-0.197 


4.33 


— 


— 


val 
c 2 


1.5 


2.10 


— 


3.29 


— 


— 


— 


„val 
C 3 


— 


5.23 


— 


1.17 


— 


19.9 


— 


„val 
c 4 


1 


— 


— 


1.5 


— 


— 


— 


„val 
C 5 


— 


2.667 


— 


— 


— 


— 


— 


sea 
c l 


— 


— 


1 


4.54 


8.19 


8.05 


— 


„sea 
c 2 


— 


1.54 


— 


1.29 


— 


— 


— 


„sea 
C 3 




2.667 












r g 
c l 




4 




4.76 


15.2 




29.3 


c 2 




2.03 




2.44 








c 3 




1.333 













Table 1: Coefficients of the state distributions parametrized in the form 
d = (a + bs + cs 2 )/(l + ds + es 2 + fs 3 + gs 4 ). 



xP g ™(x,s)^~h(s) h{z 9 ) 



^r a ;Z^s)^-h(s) S -I 2 (z s ) , (24) 
z i y 



where h(s) = exp [(—297 — 2n/)s/54], z v = 2y / 4sy/3, z s = z g = 2y/6ys, and y = hx(l/x). 
Since I„(z) — ► exp(z)/y/2nz for z — >• oo, the gluon and sea distributions grow as 
exp 2^6sha.(l/x) at small x. However, for realistic x and Q 2 values, the distributions are 
far from being asymptotic, and the power-like parametrizations (|2l|) are better suited. 

The effects of quark masses have not been taken into account in the evolution. For 
charm (and bottom), the sea distributions should therefore be modified. Different levels 
of sophistication can be used. The simplest is to assume that branchings g — > cc are 
forbidden at scales below the charm mass, while m c can be neglected in the evolution 
above threshold. In this context, m c can be thought of as an effective parameter, to 



be adjusted for taking into account threshold effects |24| . In the region of small s, the 
number of gluons increases like s, that of light sea quarks like s 2 , and that of charm sea 
quarks like s 2 — s 2 , with s c = s(max(/c 2 , m 2 ), k 2 ). For simplicity, one may thus assume 



HZi^ s) = 0{ 8 - s c ) 1 - U fT£{x, s) , (25) 




where 8(x) is the ordinary step function. 

Also original splittings 7 — > cc (and 7 — > bb) should be suppressed by mass thresholds. 
Again assuming a simple step threshold, one obtains f2' cc (x, Q 2 ; k 2 ) = for Q < m c , while 
the distribution for Q > m c is given by the substitution 

s = s(Q 2 } k 2 ) -> s(Q 2 , max(k 2 } m 2 c )) . (26) 
3.2 The anomalous distributions 

It would have been very convenient, had it been possible, to generate the full distribu- 
tions in a simple way from the state ones above. Unfortunately, the expressions are too 
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complicated to allow the dk 2 /k 2 integral in eq. (^j) to be carried out. The parametrization 
of fy> anom is therefore done separately from the ones given above. 

The /J' anom (x, Q 2 ; Qg) distributions do depend both on the Q/A and the Q /A ratio, 
and not just on a single variable. However, it is possible to separate the two depen- 
dences and to obtain a parametrization that is analytic in both Q and Qq. Hence, our 
parametrization can easily be used for a virtual photon target. 

To this end we introduce the functions 7^' qq defined by: 

7T s ) = r b ds 6Xp( 7 fe) s) . (27) 

Jo 1 — exp (— bs J 

Here s is defined analogously to s in eq. flT9|), i.e. s = (1/6) ln[ln(Q 2 /A 2 )/ln(Qo/A 2 )]. 
The /a' qq distributions can be seen as a weighted mean of the /2' qq ones; in particular 
7^' qq has a unit momentum sum, 

£ f 1 dxxTa Cm (x,s ) = l. (28) 

The information on the fraction of photons that are split by 7 — > qq is thus only carried 
by the prefactor proportional to \ii(Q 2 /Qq) so that: 

x/J— (x, Q 2 ; Ql) = ^ £ 2eJ In (g) zfTV, s ) . (29) 

The 7^' qq have been parametrized in the same spirit as the /J' q ^ ones, with due respect 
to the limiting behaviours for s — > 0. For the valence and sea parametrizations we take 
the same functional form as in eq. fl2~T|), but for the gluon we choose 

xT^ix, s ) = cf x c " (1 - x 2 ) c i xf^f(x) . (30) 

This form is not exactly like that for the gluon state-distribution, rather ((1— x)(— lnx))° 3 
is replaced by (1 — x 2 )° 3 . We have attempted to use the same form in both cases, but 
then failed to obtain reasonable fits. The limit s — > is again exact. The coefficients Cj 
are collected in table 0. 

The small-x behaviour is again better described by a powerlike distribution over the 
relevant x ranges, even though asymptotically the growth is weaker, for example 

x tf—<*. Q\Ql) « £ (E 2e?) in g ^ exp (-^^ .) /, W . (31) 

The charm (and bottom) seas are suppressed; close to threshold by a relative amount 
1 — (s c /s) 3 , where s c = s(max(Qo, m l), Ql)- (Note that this factor approaches unity faster 
than the factor in eq. ( p5|) , as a consequence of having a spectrum of initial branchings at 
different k 2 scales.) For Q > m c , the charm sea is therefore assumed to be given by the 
normal sea multiplied by this factor. 

From the above formulae, the full distribution is built up by summing the contributions 
from the allowed flavours as already specified in eq. (p9|). Branchings 7 — > cc (and 7 — > bb) 
do not contribute below threshold. Above threshold, the charm contribution is obtained 
by substituting max(Qg,m 2 ) for Qq in the logarithmic prefactor and in the definition of 
s. 
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a 


b 


c 




e 


/ 


9 


«val 
c l 


1.5 


2.49 


26.9 


— 


32.3 


— 


— 


val 
c 2 


1.5 


-0.49 


7.83 


— 


7.68 


— 


— 


„val 
C 3 


— 


1.5 


— 


-3.2 


7. 


— 


— 


„val 
c 4 


1 


— 


— 


0.58 


— 


— 


— 


„val 
C 5 


— 


2.5 


— 


10. 


— 


— 


— 


sea 
c l 


— 


— 


0.333 


4.90 


4.69 


21.4 


— 


„sea 
c 2 


— 


-1.18 


— 


1.22 


— 


— 


— 


~sea 
c 3 




1.22 












r g 
c l 




2 




4. 


7. 






r g 
c 2 




-1.67 




2. 








c 3 




1.2 













Table 2: Coefficients of the full distributions parametrized as in table [371 . 



3.3 The heavy-flavour contribution 

The above treatment of heavy-flavour branchings 7 — > cc, bb is not unreasonable for jet 
production in 77 collisions of real photons, Q 2 = P 2 = (with the role of hard scale here 
taken by the transverse mass of the quarks). However, in the deep- inelastic-scattering 
region, Q 2 ^> P 2 « 0, the kinematics constraint W 2 = Q 2 (l — x)/x > <im 2 implies that 
the distribution can only be non- vanishing for x < Q 2 /{Q 2 + 4m^). The large charm- 
quark charge ensures that the charm contribution is very significant for the shape of the 
F1{ x tQ 2 ) distribution. 

For the charm contribution to F2, we have therefore based ourselves on the leading- 
order "Bethe-Heitler" cross section for 7*7* — > cc. (An additional charm contribution 
comes from our sea parametrization, i.e. from g — > cc branchings.) The full expression 



is given in f25fl , but for our applications it is sufficient to use the approximation in |2 
valid in the limit Ax 2 P 2 < Q 2 : 

xc(x, Q 2 ) = ^3 (^) 2 x {ft [6x(l -x)-l] 



+ log 



PcVc 



1 - Pen 

2x 2{3 cVc 
Q 2 l-(3 2 cV 2 



x 2 + (1 — x) 2 + r c x(l — 3a;) — -r 2 x 2 
(2 - r c )m 2 c x - P 2 x] \ , (32) 



with 



Ami „ / 4m? / Ax 2 P 2 



c 



c Q2 , & = yi-^, ^ = v 1_ "^~' (33) 

and correspondingly for b. 

Some distance above the threshold, the Bethe-Heither formula gives essentially the 
same contribution to F 2 as does the standard formulae presented in the previous section. 

In the calculations, we have used fairly small c and b masses, m c = 1.3 GeV and 
mb = 4.6 GeV, so that the charm contribution is large. On the other hand, the J/ip and 
T states are not included in the VMD sum over vector mesons. To first approximation, 
these two effects should therefore cancel. The approach can be motivated by standard 
duality arguments. 
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4 The non-perturbative part 



4.1 Normalization 

We now turn to the hadronic input distributions fa' V ( x ) — fa' V ( x > Qoi Qo)j °f- ec l- (0)- 
Usually these are estimated by the pion PDFs using VMD and the additive quark model. 
Of course, the PDFs of the qq "bound states" of the photon need not be the same as 
those of real vector mesons. Moreover, the PDFs of the short-lived //'-meson may well 
differ in shape from those of the long-lived pion. Therefore, in contrast to the GRV 0] 
and AFG [K| parametrizations, we do not approximate the hadronic input distributions 
fa ,V ( x ) by the pion input distributions fa( x ), but rather constrain the shape of the input 
distributions from data. 

We do, however, use VMD to fix the normalization of the input distributions. Indeed, 
for consistency with eq. ([!]) we cannot allow for an extra "K-factor" as is introduced in 
refs. [|, |10J if the input scale Qo corresponds to the simple p°, u>, VMD model. Based 



on eqs. (||), (||) and ( |TTD one has 

r<(x,Q 2 )-r<^(x,Q 2 ) = 

E ^fl y (*, q\ 05) + ^ E K C % f«^ x > ^ ^ • ( 34 ) 

y Jy Z7r q J Q K 

The analogous decomposition ([]]) of the 7p total cross section tells us that Qo should be 
of the order of 600 MeV if the sum over vector mesons includes only the three lowest-lying 
states V\ = p°, a), and 0. On the other hand, for larger values of Qo more vector mesons 
V n have to be included: 



E f \v) E E -^\v n ) . (35) 

V=p°,ui,<j) J V n=l V=p°,Lu,<f> J Vn 

Since the VMD couplings of such higher-mass states are poorly known, mass effects be- 
come more important. Moreover, the input distributions will anyhow be fitted, so we may 
approximate 

N( - Qo) Ana „ 47ra ~ „ 

E E E ^Pa v ^)- (36) 

n=l V=p°,uj,4> JV„ V=p°,(j,4> JV 

In principle, for a given Q , K = K(Q ) could be determined by repeating the jp total 
cross-section analysis. In the following we will consider two cases, (i) Qo — 0.6 GeV and 
K = 1 fixed from the jp analysis, and (ii) Qo = 2 GeV, where K will be fitted together 
with the input distributions f(x) to the F% (x, Q 2 ) data. 

The normalization of the PDFs discussed so far still leaves us freedom of how to add 
the vector mesons. For example, in the case of simple VMD: 



|7)VMD 



e 2 e 2 / 9 9\-V2 / , . , e 



E t ^ = \ 75 + « ( e * + e *)~ ( e - l uTI > + ^ l dd )) + T ^ ■ ^ 

V=p,uj,4> J y \ Jp Ju J<t> 



The coherent superposition of vector mesons is obtained for e u = 2/3, = —1/3. 
A completely S^-symmetric coherent superposition would give f p /f u = 1/3 and 
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a 


Incoherent 


Coherent 


Coherent (naive S^) 


u-val 


(£ + <£)/2 ^0.248 


4(^ + ^)/5 ^0.397 


8^/9 « 0.404 


u-sea 


r> T T 1 

9i + 9i + 9i -0.551 


n — t n 

9i + 9t> + 4 ~ 0.551 


4^/3 « 0.606 


d-val 


(g 2 + gl)/2 ^0.248 


(^ + ^)/5 -0.099 


2^/9 « 0.101 


d-sea 


9i + 9t + 9i ~0.551 


^ + ^ + ^ ^0.551 


4^/3 « 0.606 


s-val 


^ « 0.054 


^ « 0.054 


2^/9 « 0.101 


s-sea 


gl + gl + gl - 0.551 


gl + gl + gl « 0.551 


4^/3 « 0.606 


gluon 


£ + 9% + 9* « 0.551 


+ £ + -0.551 


4^/3 « 0.606 



Table 3: The coefficients c a of (|38|) for coherent and incoherent superpositions of vector 
mesons. Here g v = 4tt/ f v . 



fp/f<P = —V2/3. This is close to the experimental numbers, but in particular (f> pro- 
duction seems suppressed compared to the above expectation. This is not unreasonable 
in view of the larger s quark and <fi meson masses. We therefore keep the VMD couplings 
at their measured values, which we also use for the total cross section in eq. ([!]). An 
incoherent superposition of u and d quarks corresponds to e u = e& = 1. At long time 
scales, e.g. in "elastic" processes such as 7p — > Vp, the p° and u vector mesons contain 
equal amounts of uu and dd, i.e. the coherence at the 7qq vertex is broken. It is therefore 
conceivable to have either physics scenario, or some intermediate thereof, and both have 
been used in the literature [1|, H H3 • Here we favour the coherent superposition of u and 
d quarks, in line with the argument that hard processes probe short time scales. 

For the input distributions we assume an S'LVsymmetric sea distribution s(x) and 
denote by v(x) the input valence distribution such that £ = J2q( X( l( x ) + xq(x)) = 
2v(x) + 6s(x). The coefficients of the various hadronic input distributions /a' VMD (x) 
are summarized in table |3| for the three possibilities of photon decomposition discussed 
above, and for K = 1. We have used an obvious notation, e.g. 

/Z/vT D (^ Qo, Qo) = «em C a , va l v(x) , (38) 

and our preferred choice corresponds to the middle column in table For completeness, 
we list also the expression of F 2 in terms of s(x) and v(x) for this choice 



a: 



f^' vmd (x, Q 2 ) 



cm 



34 Utt 4tt\ 2 4vr 

45 [tf+J! yyJI 



6 v=p,u),<t> Jv 



= 0.3875 x v(x) + 0.735 x s(x) . (39) 
We choose the following ansatz for the input valence distribution 



xv 



(x) = N v x a " (1 - x) bv 1 + C v y/x (40) 



and analogous forms for the sea and gluon distributions. Note that, since we fix the nor- 
malization of the hadronic input distributions using VMD, we can exploit two constraints, 
namely the number of valence quarks and the momentum sum rule 

[ 1 dx2v(x) = 2 ; dx V x Ux) = 1 . (41) 
Jo Jo 
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4.2 Direct contribution 



In contrast to the proton structure function F 2 , the photon structure function F 2 receives 
a direct contribution beyond LO accuracy, symbolically 



Cq + /g 



a + c 1 . 



(42) 



(Note that summation over q does not automatically include antiquarks as well; the 
additional q contribution is reflected in a factor 2 in many of the subsequent formulae.) 
In fact, at the relatively large x values currently probed, the inclusion of the direct term 



C 7 (x) 



o 



2n 



3E 



4 



x 2 + (l-x) 2 



In 



x 



X 



+ 8x(l 



1 



(43) 



is much more important than the modification in the Q 2 evolution when going from 
LO to NLO. Inclusion of C 7 defines the MS scheme, while C 7 = defines the DIS 
scheme. At any order in perturbation theory, the theoretical prediction for F 2 is scheme- 
dependent. This factorization-scheme dependence is reduced when including higher-order 
contributions (i.e. when going from LO to NLO to NNLO etc.). Many cross sections in 
7P and 77 are still calculated in LO accuracy. In order to investigate the factorization- 
scheme dependence associated with a LO calculation we shall fit PDFs in LO for both 
schemes 

DIS: F2(x,Q 2 ) = 2j2elf<(x,Q 2 ) 

1 

MS: F 7 (x,Q 2 ) = 2^ e 2/ 7 (a;,Q 2 )+C 7 (x) , (44) 
q 

although formally C 7 is of NLO. However, we only take that part of C 7 that is universal 



Til , namely 

C^{x) 



2tt 



3 2 4 X 



X 



X) 



In — h 6x(l 

x 



x 



(45) 



This modified expression follows from a careful analysis of the ladder diagrams, but can 
also be obtained by simply imposing a t (and u) cut t on the box diagram to separate 
the perturbative (point-like) from the non-perturbative (hadronic) parts resulting in 



^7, box 



O 



2n 



x,Q 2 ) 

-3E 



e 4 x 



x 



X 



In 



to 



- + [6x(l - x) - 1] ( 1 



2t ( 



(46) 



with t r 



Q 2 /x. Upon dropping the higher-twist terms (oc l/Q 2 ), eq. 



]) yields 

eq. fl45D . There is no point in keeping the higher- twist contribution in C 7 when neglecting 
the further (unknown) higher-twist contributions to F 2 . Since charm (and bottom) are 
included by the Bethe-Heitler formulae, they are not affected by the considerations in this 
section. Hence the quark sum here only runs over u, d and s quarks. 



13 



4.3 Low-Qo fit 

The F~2 data are not yet precise enough for the QCD scale parameter A to be fitted. 
Therefore we fix (one-loop) A(n/ = 3) = 230 MeV corresponding to A(nf = 4) = 200 MeV. 
Moreover, for K = 1, Q = 0.6 GeV is fixed from our 7p analysis. In the DIS scheme, 
the simplest choice, a valence-like input (N s = and g(x) oc v(x) with v(x = 0) = 0), 
results in a reasonable x 2 only for a considerably smaller Qo ~ 0.4 GeV. Therefore we 
allow for sea-like input distributions. As the data do not require non-zero values for the 
parameters a s and Cj, cf. eq. fl4"0D, we simply put these parameters equal to zero. The b g 
and b s are poorly constrained, and have therefore been assumed to be given by b g = b v + l, 
b s = b v + 3. Finally, a g = a v /2 is not inconsistent with the data and has been fixed. Then 
we find in the DIS scheme 

set SaS ID input : 

DIS: Q = 0.6GeV , A = 0.2 GeV 
x v(x) = 1.294x a80 (l - xf 76 
x s(x) = 0.100(1 -x) 3 ' 76 

x (7(x) = 1.273x ' 40 (l-x) L76 (47) 

at a x 2 °f 141 for 71 i 7 ^ data points. Note that we have included also the low-Q 2 points 
(Q 2 = 0.71 GeV 2 ) in the fit. At the input scale Q = Qo, the photon momentum is split 
according to about 5:1:2 between the (two) valence quarks, the (six) sea quarks, and 
the gluon. 

In the MS scheme the fit favours a small sea and we find 

set SaS 1M input : 

MS: Q = 0.6GeV , A = 0.2 GeV 
x v(x) = 0.8477x a51 (l - x) 137 
x s(x) = 

x g(x) = 3.42x a255 (l - xf 37 (48) 

at a x 2 of 136, where the photon momentum is now split according to about 7 : : 13 at 
the input scale Q = Qq. Both Q 2 -evolved PDFs are parametrized as 

xf2' YMB {x,s) = c?x c 2 (l-x) c 3 [-lnx] c 4 + c£x c 6 (l-x) c 7 a = q, val, q, sea, g , (49) 

and the coefficients Cj are given in tables |4] and [5], respectively. The coefficients C5-C7 are 
only listed in those cases where the second term was included. 

Figure ^ shows F^x, Q 2 ) data compared with the predictions of our sets. As explained 
above, these consist of a VMD part (different for the two sets), an anomalous part for 
u, d and s quarks (in common), Bethe-Heitler terms for c and b production (also in 
common), and, for set 1M only, the C 1 term for u, d and s quarks. This subdivision 
is illustrated in Fig. ^|, for one specific Q 2 scale. With the rather low value of Qo, the 
anomalous contribution quickly becomes the dominant one. At higher Q 2 and smaller x, 
also the charm contribution is important. In the MS scheme the hadronic distributions 
vanish faster with x than in the DIS scheme, both since the C 7 term is negative at large 
x and since the (1 — x)-power is larger. 

A reasonable description is obtained for all Q 2 , although there might be a slight 
incompatibility between the low- and high-Q 2 data, which is also reflected in the not quite 
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a 


b 


c 


d 


e 


/ 


9 


K 


c i 


1.294 








0.252 


3.079 






-iva.1 
c 2 


0.80 


-0.13 














^val 
c 3 


0.76 


0.667 














„val 

c 4 




2.0 














c l 


0.1 




-0.397 


1.121 




5.61 


5.26 




C 2 






-7.32 






10.3 






C 3 


3.76 


15.0 


12.0 




4.0 








c 4 


















cf 




7.90 






5.50 






5.16 






—1.9 






3.6 








r g 

c 3 


1.3 
















r g 
C 4 


0.5 


3.0 














r s 

c 5 


1.273 














10.0 


r g 
c 6 


0.4 
















r g 
c 7 


1.76 


3.0 















Table 4: Coefficients of the hadronic distributions evolved from (|47|) , parametrized in the 
form fl49|), where q = exp(— Ks)(a + bs + cs 2 + ds 3 )/(l + es + fs 2 + gs 3 ). 





a 


b 


c 


d 


e 


/ 


9 


K 


„val 
c l 


0.8477 








1.37 


2.18 


3.73 




val 
c 2 


0.51 


0.21 














„val 
C 3 


1.37 
















„val 
c 4 




2.667 














„sea 
c l 




0.842 






21.3 


-33.2 


229. 




„sea 
c 2 


0.13 


-2.90 






5.44 








sea 
C 3 


3.45 


0.5 














„sea 
c 4 


2.8 
















cf 




24.0 






9.6 


0.92 


14.34 


5.94 


4 


-0.013 


-1.8 






3.14 








c 3 


2.37 


0.4 














r S 
c 4 


0.32 


3.6 














r s 

c 5 


3.42 














12.0 


c 6 


0.255 
















r g 


2.37 


3.0 















Table 5: Coefficients of the hadronic distributions evolved from (|4*8"|), parametrized in the 
form fl49|), where C; = exp(— ns)(a + bs + cs 2 + ds 3 )/(l + es + fs 2 + gs 3 ). 



15 



optimal x 2 - This could indicate that higher- twist contributions are not negligible at low 
Q 2 . One can imagine a situation in which higher- twist contributions significantly affect 
F 2 , but where the leading-twist evolution of PDFs is still valid down to Qq = 0.6 GeV. In 
fact, including only the F^x, Q 2 ) data above Q 2 = 4 GeV 2 improves the fit to a x 2 of 60 
(78) in the MS (DIS) scheme at now 55 data points. However, in both cases a very large 
sea component is needed. Therefore we do not provide such a parametrization but rather 
one where we also take a large Qq = 2 GeV. 

The distributions ([IT]) vanish do truly hadronic distributions. In order to 

check whether the data prefer a harder (point-like) component we added a term N v d v x to 
eq. ( fHf ) but found no improvement in x 2 ■ We checked also the normalization of the VMD 
couplings by allowing the values g 2 in table |3| to be multiplied by an overall constant K. 
Interestingly enough, the data do not really require values for the VMD couplings that 
deviate notably from their values used in the description of 7p cross sections. The fit gives 
K — 1.17 with only a slightly better x 2 (136 compared to 141). This 17% deviation of K 
from unity is within the range of uncertainty of the gf. For example, using / 2 /47r = 2.02, 
as extracted from the leptonic width alone, gives a 9% increase compared to the value we 
use, namely 2.20, which is the geometrical mean between this "leptonic" value and the 
value extracted from 7p — > p°p. 



4.4 High-Qo fit 

Alternatively we fit input distributions starting at a "typical" deep-inelastic Qq value, 
namely Qq = 2 GeV. Fixing again A 4 at 200 MeV, as well as a s = a g = Ci = 0, b s = 4, 
b g = 2, we find in the DIS scheme K = 2.422 (corresponding to additionally allowed 
vector mesons between 0.6 GeV and 2 GeV) and 

set SaS 2D input : 

DIS : Q = 2 GeV , A = 0.2 GeV , K = 2.422 
Kx v(x) = 1.00 [a; - 46 (l - x) 064 + 0.76a;" 
Kx s{x) = 0.242(1 - xf 

Kx g(x) = 1.925(1 - x) 2 (50) 

with x 2 — 59 for 55 F 2 data points above Q 2 = 4 GeV 2 . In this case the photon momentum 
is split as about 5:1:2 between valence quarks, sea quarks, and gluons, where the 
"hadronic" term and the "hard" term in eq. ( |50"D contribute each about 50% to the 
valence momentum. As can be seen from Fig. |l|, this fit describes the high-Q 2 data better 
than the fit starting at Qo = 0.6 GeV, however, at the expense of not accommodating the 
\ow-Q 2 data. 

As in the case of the low-Q 2 fits, the fit in the MS scheme again favours softer distri- 
bution functions and we find 



set SaS 2M 
MS : 
Kx v(x) 
Kx s(x) 
Kx g(x) 



input : 

Qo = 2 GeV 



A = 0.2 GeV 



K = 2.094 



x°- 50 (l 



1.168 
0.209(1 -x) 4 
1.808(1 -x) 2 



x 



,2.60 



+ 0.826x 



(51) 
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a 


b 


c 


d 


e 


/ 


9 


K 


n vai 
c l 


1.0 


0.186 


— 


— 


-0.209 


1.495 


— 


— 


val 
c 2 


0.46 


0.25 


— 


— 


— 


— 


— 


— 


„val 
C 3 


0.64 


0.14 


5.0 


— 


1.0 


— 


— 


— 


„val 
C 4 


— 


1.9 


— 


— 


— 


— 


— 


— 


„val 
C 5 


0.76 


0.4 


— 


— 


— 


— 


— 


— 


val 
C 6 


1.0 


— 


— 


— 


— 


— 


— 


— 


„val 
C 7 


— 


2.667 


— 


— 


— 


— 


— 


— 


cf a 


0.242 


-0.252 


1.19 


— 


-0.607 


21.95 


— 


— 


~sea 
c 2 


— 


— 


-12.1 


— 


2.62 


16.7 


— 


— 


~sea 
c 3 


4.0 


— 


— 


— 


— 


— 


— 


— 


~sea 
c 4 


— 


1.0 


— 


— 


— 


— 


— 


— 


cf 


1.925 


5.55 


147. 




-3.59 


3.32 




18.67 


c! 




-5.81 


-5.34 




29.0 


-4.26 








2.0 


-5.9 






1.7 








cf 




9.3 






1.7 









Table 6: Coefficients of the hadronic distributions evolved from (|50|) , parametrized in the 
form (|4"9"|), where c« = exp(— ks)(q + bs + cs 2 + ds 3 )/(l + es + fs 2 + gs 3 ). 



with the same \ 2 an d a slightly larger gluon- momentum fraction (29% compared to 27%). 
The comparison with data and the low-Q 2 fit is given in Fig. [I]. 

The breakdown of the two sets component by component is again shown in Fig. |^. 
Owing to the larger Qq scale, the anomalous component is small at the lower Q 2 range. 
Note that the Bethe-Heitler contribution to cc production remains unchanged, since we 
do not include an increased amount (or, indeed, any amount at all) of J/ip production in 
the VMD component and therefore are allowed to maintain the same m c = 1.3 GeV as 
was used for the low-Qo fits. 

The distributions have again been parametrized according to eq. ( ^9|) , with coefficients 
given in tables |5] and [7], respectively. 

The current uncertainties in the PDFs of the photon are displayed in Fig. |3], which 
shows the u-quark and gluon PDFs. For the quark distribution, the high-Q 2 distributions 
SaS 2D and 2M are larger than the low-Q 2 ones ID and 1M, simply reflecting the corre- 
sponding hierarchy in K7. At small x the opposite behaviour holds true since the small-x 
rise increases with the length of the Q 2 evolution. Hence the low-Q 2 parametrizations 
also yield larger gluon distributions at small x. 

The effect of a non-zero target virtuality is displayed in Fig. [|. Increasing P 2 not only 
lowers the normalization but also changes the shape of the distribution. As expected, the 
gluon distribution vanishes faster than the quark distribution as P 2 — > Q 2 . 

4.5 The TPC/27 parametrization and comparison with FKP 



Experiments often use the TPC/27 parametrization of i<7 (x, Qq) determined at Q 
0.71 GeV 2 to estimate the hadronic component of F^{x, Q 2 ) 



1 i^ VMD (x, Q 2 ) « — Fj(x, Q 2 )|tpc/ 27 ( 52) 
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a 


b 


c 


d 


e 




9 


K 


«val 
c l 


1.168 


1.771 


29.35 


— 


— 


— 


— 


5.776 
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Table 7: Coefficients of the hadronic distributions evolved from (|5ll) , parametrized in the 
form (091), where q = exp(— ks)(q + bs + cs 2 + ds 3 )/(l + es + fs 2 + gs 3 ). 



= (0.22 ± 0.01) a; a31±a02 (1 - a;) ' 95 + (0.06 ± 0.01) (1 - x) 2 - 5±1A 

The full F^(x, Q 2 ) distribution is then obtained by adding the FKP parametrization JZIJ 
of the anomalous part of to eq. (0). The only free parameter^ in this aproach is the 
cutoff scale pto separating the perturbative part of from the non-perturbative one. 
Moreover, pro is often varied in the figures of F 2 7 (x,Q 2 ) or f°fdxF2(XiQ 2 ) to illustrate 
the strong dependence on this arbitrary scale, implying little or no sensitivity of the data 
to the QCD scale A. This approach to compose from a hadronic and a perturbative 
part can be criticized on three counts: 

1. The hadronic part of i 7 ^, eq. (|52"[), is not evolved in Q 2 . 

2. The strong correlation between the scale pro and the size and shape of the hadronic 
part is neglected. 

3. The FKP formula is a poor parametrization of the anomalous part. 
Let us discuss these points in more detail. 

The first point is obvious. In order to investigate how well the F% data can be described 
with eq. (|52"D, one has to consider this equation as the non-perturbative input to F^(x, Q 2 ) 
at the scale Qq = 0.71 GeV 2 , whereafter it evolves with Q 2 . Note that, for the Q 2 
evolution, the input distributions have to be known separately for the valence, sea, and 
gluon distributions. It is sensible to assume that these hadronic input distributions can be 
interpreted within generalized vector-meson dominance (GVMD), in which case we can 
exploit eq. (0). Associating the first term of eq. ([52]) with valence quarks, we need a K- 
factor of K = 1.411 compared to simple VMD, eq. (p9|). This represents the contribution 
of additional vector mesons that have to be included in eqs. (|ll|) and ( p9|) because of the 
larger continuum cutoff Qq = 0.84 GeV compared to the ~ 0.6 GeV preferred for simple 
VMD. 



2 The values of the light-quark masses are needed as well, but these are supposed to be known. 
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The only free parameters are then A, a g , and b g . 
scheme yields A = 0.116 GeV and 



The fit to the F 2 data in the DIS 



Kx v(x) 
Kx s(x) 
Kx g(x) 



„0.31/ 



\0.95 



0.082(1 - x) 2 - 5 
2.586s a31 (l-x) - 95 



(53) 



The photon momentum is approximately split as 3 : 1 : 6 between valence quarks, sea 
quarks, and gluons at the input scale. In Fig. || we show the hadronic part of F 2 at Q 2 = 4 
and 100 GeV 2 as well as at the input scale 0.71 GeV 2 . The effect of the Q 2 evolution is 
non-negligible: a clearly visible softening at large x and an increase at small x is observed. 
The fit is, however, a rather bad representation of the data, in particular at high Q 2 . This 
is reflected in a high x 2 value, \ 2 = 181. We therefore do not provide a parametrization 
for the distributions Q 2 -evolved from eq. fl52]). 

Concerning the second point, we recall that, by definition, the perturbative part of 
F 2 (x, Q 2 ) vanishes at the Qo (or pto) scale separating it from the non-perturbative part 
of F 2 . Hence, any change in this scale has to be accompanied by a corresponding change 
in the size and shape of the hadronic input distributions at this scale. The correlation is 
obvious: the larger the Qo, the smaller is the range of perturbative evolution, and hence 
the bigger is the non-perturbative input needed. If higher-twist effects can be neglected, 
the total photon momentum carried by partons should be independent of Qo- Then 
eqs. (^) and ( fT2p yield the following dependence on Qo for the normalization K = K(Q ) 
defined in eq. (|36|): 



K(Q' ) = K(Q ) + 



T e 2 



V=p,u>,- 



Ql 

MUV^ Ql 



in 



K(Q ) + 0.770 In ^ . 

Qo 



(54) 



Also a correlation between the separation scale and the x-shape is to be expected, 
with harder hadronic input distributions required for larger Q Q scales. This could be ar- 
gued on the basis of GVMD, or simply by the fact that parts of the (hard) perturbative 
7 — > qq splittings have to be included in the "hadronic" input distributions for larger Q . 
In any case, the separation scale is an arbitrary parameter and hence the result for F 2 
must not depend on its value, provided it is changed within a range in which leading-twist 
perturbative QCD describes the physics. As the comparison of our various parametriza- 
tions demonstrates, even the variation of the scale between 0.6 GeV and 2.0 GeV has not 
a major influence on the final F 2 result. (The differences that are there come from them 
being fits to different data sets, i.e. data with Q < 2 GeV included or not. One can do 
a Qo = 0.6 GeV fit to Q > 2 GeV data only, and then obtain close agreement with a 
Qo = 2 GeV fit.) 

As to point 3 above, a solution of the perturbatively calculable part of the quark 



distribution functions of the photon has been obtained by FKP [20]: 



(x,Q 2 



a. 



2 x 2 + (l-x) 2 
2tt q x c + Cf(x) 



T 



1-4^ 



To 
T 



1+Cf(x)- 



(55) 



where 



C 



8 

33 - 2N 



f 



/(*) 



2 In 



x 



x 



2 X 
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T 



i ^max 

111 A 2 


^max 


= 91 

X 


- Ill —T 

A 2 


to = 


2 

m T0 _ 
1 — X 



VP? + Pyg + X(l — X)P 2 



X 



(56) 



Here to is the cutoff in the integration over the virtuality of the t-channel and w-channel 
quark propagators and P 2 the target mass (virtuality of the probed photon). In Fig. ||a 
we compare our anomalous u-quark distribution function at Q 2 = 4GeV 2 , P 2 = 0, and 
Qo — 0.36 GeV 2 with eq. fl55|), where we have taken the FKP value for the quark mass 
(m q = 0.3 GeV) but adjusted pro such that m\ Q = Qq. 

It might seem surprising that the two distributions do not at all agree with each 
other. The difference arises mainly from the inclusion of the 1/x and 1/(1 — x) factors 
in eq. ( ]5^| ) (1/x in t max and T, 1/(1 — x) in t and T ). These terms are of kinematic 
origin and specific to the quark-parton-model result (calculation of the box diagram). 
However, their inclusion in the leading-twist QCD evolution equations lacks justification 
in perturbative QCD: the leading-logarithmic approximation sums only the logarithms in 
Q 2 /A 2 and makes no statement about the summation of lnx or ln(l — x) terms (beyond, 
of course, the double-leading-logarithmic approximation). In other words: although the 
factorization scale is arbitrary and may well differ from its "natural" value Q 2 , it must 
not depend on x in order not to spoil factorization. 

On the other hand, it is satisfying to observe that the FKP parametrization agrees 
almost perfectly with ours at large x if we take t max = Q 2 and t Q = m^ Q in eq. (|56|). 
Differences at small x arise because eq. (^) was derived in the valence approximation, 
hence xf q given in eq. Q5"5] ) vanishes as x — > while the full solution gives xf q rising at 
small x. 

Of course, one could argue that the FKP formula is not meant as a parametrization 
of quark distributions but as part of an F^x, Q 2 ) parametrization, in combination with 
VMD and direct terms. (After all, it is the FKP F2 expression that is being used in 
the experimental analyses.) In this context, In a; and ln(l — x) terms may well be redis- 
tributed from the direct term into the "quark distribution" if this provides a convenient 
approximation to the full solution. The FKP prescription for F^ is 



Fr n ° m (x, Q 2 ) = 2 £ e\x \ K"(x, Q 2 ) 

2xm 2 — x(l — x)P 2 



q 



In 



6x(l — x) — 1 + 



q 



2 

m T0 



(57) 



Figure [f]b compares the u-quark component of eq. ([57]) with our parametrization of the 
anomalous u-quark distribution function, to which we add the u-quark component of the 
universal direct term eq. (H3), i.e. 



^fim . / 2 



C2(x) = -^3 - ) x 



2tt V3 



2 



x 2 + (1 - x) 2 } In - + 6x(l - x) - l) . (51 
J x J 



The observable discrepancy is, in fact, not unexpected as no ln(l — x) term is included in 
eq. flSSp. However, a discrepancy remains even when the ln(l — x) factor is removed from 
the definition of T in the FKP formula. 

The reason for this is the extra m 2 /my term in eq. (|57|) . Indeed, when we take m q = 
and t = Qq in eqs. (|56|) and (|57D, we find agreement at large x between our anomalous- 
plus-direct u-quark distribution function and the FKP one, Fig. ||c. Alternatively we 
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also expect agreement for the comparison of the standard FKP result, i.e. where t = 
Qq/(1 — x), if we change from the expression (f)8|) for the direct term to the standard 
direct term, eq. 



C&x) = ^3n i f[i 2 + (l- x) 2 ] In - — - + 8x(l - x) - l) , (59) 
Z7T \3J L L J X J 

which includes the ln(l — x) term. This is, however, not the case. The difference can 
be traced back to the difference in the coefficient of the x(l — x) term in eq. ([57]) and in 
eq. ([591). Indeed, we obtain agreement with the standard FKP result at large x if we add 
to our parametrization of the u-quark distribution function the following direct term: 



ftftm ~ / 2 



C2(x) = — 3 



2tt V3 



) x | [x 2 + (1 - x) 2 ] In ^— - + 6x(l - z) - 1 j . (60) 



In conclusion, putting m q equal to zero and changing the 6x(l — x) term into 8x(l — x) 
in eq. (|57D , the FKP parametrization describes rather well the leading- twist, leading-order 
anomalous part of F 2 (x, Q 2 ) at large x in the (standard) DIS scheme. However, a number 
of problems remain: 

• The unmodified FKP expression (^7j) includes extra terms that have no correspon- 
dence in perturbative QCD: they are neither leading- twist NLO terms nor higher- 
twist contributions. These terms depend on an additional, unphysical parameter, 
the light-quark masses. 

• The FKP result is obtained in the valence approximation, and thus does not re- 
produce the correct small-x behaviour of F 2 (which, actually, matters already for 
x < 0.5). 

• A hadronic part of F 2 matched to the FKP parametrization is not available. The 
recipe "TPC/27-plus-FKP" is inconsistent: since the FKP parametrization needs 
a Pro < 0.5 GeV to describe the high-Q 2 data, FKP is already large at Q 2 = 
0.71 GeV 2 , where the TPC parametrization alone is supposed to describe the data. 

• Since the FKP parametrization in eq. (^) mixes kinematic terms between the 
evolved quark distribution functions and the direct term, it is not a parametrization 
of the anomalous quark distribution functions of the photon (at least not with de- 
fault parameters). The quark distribution functions do not show the correct small-x 
behaviour and a parametrization of the gluon density is missing. 

We have attempted to provide an alternative where these problems are resolved. Our 
quark and gluon distribution functions obey the leading-logarithmic evolution equations 
down to very small x (~ 10~ 4 ). The only free parameter besides A is the scale Qo- 
Moreover, for two representative values of Qq (0.6 and 2.0 GeV), our anomalous PDFs 
are supplemented by parametrizations of the hadronic PDFs. These have been obtained 
by fitting the photon structure function F 2 (x, Q 2 ) to the available data. 



5 Concluding remarks 

In this paper we have shown how the PDFs of the photon can be decomposed into pertur- 
bative and non-perturbative components. Since this separation and the constraints on the 
non-perturbative part is the central result of our paper, we briefly repeat the main ideas 
of our approach and show its connection with an integral representation of the photon 
structure functions. Our starting point is the fact that a PDF, being the solution of an 
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inhomogeneous evolution equation, can always be written as the sum of two terms: a 
particular solution f2' PT (x, Q 2 , Qq), which vanishes at Q 2 = Q 2 ,, and a general solution 
/J ,NP (x, Q 2 , Qq) of the corresponding homogeneous evolution equation, that needs an in- 
put at Q 2 = Qq. For Qq sufficiently large, such that higher-twist (1/Q 2 ) effects become 
negligible, the anomalous distribution function defined in eq. (|8|) is a particular solution. 
Hence we can identify f2' PT (x,Q 2 ,Q 2 ) = /J' anom (x, Q 2 , Q 2 ). 

To constrain /J' NP (x, Q 2 , Qq) we first observe that the anomalous distribution function 
of a parton a within the photon is the convolution of two factors. The first one, P^(k) = 
(a em /27r)2e 2 d/c 2 /fc 2 gives the probability for the photon to branch into a qq state at some 
(perturbatively high) scale k. The second factor is the distribution function of the parton 
a within this qq state. This "state" distribution function obeys the homogeneous evolution 
equation (with the calculable input (0)). 

This factorization of /J' PT suggests a similar one for /J' NP (a;, Q 2 , Qq). Below the 
scale Qo, the 7 — > qq transition can no longer be calculated perturbatively but may be 
approximated by fluctuations of the photon into vector mesons with probabilities Pyik 2 = 
m v ) = 47ra em / f v . Hence the analogue of the fc-integral above Qo of the continuous 
spectrum of perturbative ( "anomalous" ) states is the discrete sum of vector mesons, with 
the state distribution functions f2' qq replaced by the PDFs f2' v of the qq "bound states" 
V of the photon. As long as the shape of these distribution functions are treated as free 
parameters (to be fixed from data), /J' NP = /7' VMD defined in eq. (|TT[) is indeed a general 
solution of the homogeneous evolution equations. 

Our decomposition ([34]) also follows from the representation of the (moments of the) 
photonic PDFs as a dispersion integral in the photon mass (P 2 is the photon virtuality) 

m 

rQ 2 HP 

r a (n, Q 2 , P 2 ) = jT j^TpjPain, Q 2 , k 2 ) . (61) 

Rather than describing the dispersion integral as the difference between a "point-like" 
part (contribution from the upper limit) and a "hadronic" part (contribution from the 
lower limit), it is more natural to separate short- distance and long-distance parts by a 
scale Qo since the weight function p a possesses the scaling-violation pattern typical of 




ordinary hadronic PDFs p7 |. Our solution follows from the ansatz 
Pa {n, Q 2 , k 2 ) = ^e 2 /™>, Q 2 , k 2 )Q(k 2 - Q 2 ) 

7T 

(62) 

The scale Qo is arbitrary; hence the Qo dependence of the anomalous part must be 
cancelled by that of the hadronic one. For example, in the generalized vector-meson- 
dominance ansatz ([36]) that we are using, the relation between the two parts is given by 
eq- O- 

For a given, chosen Qo, the only unknown (non-perturbative) pieces of the photon 
structure functions and the PDFs of the photon are the valence, sea, and gluon input 
distributions and the value of K in eq. (0). In order to cover the uncertainties associated 
with our approach we have presented two extreme analyses. In the first we restricted 
VMD to the well-established p°, u, <fi states only (i.e. K = 1) and fixed Q = 0.6 GeV 
as obtained from our 7p analysis. Here the main theoretical error arises from the use of 
perturbation theory down to rather low values of Q 2 . The spirit of the second analysis 
was opposite: take Qo well within the perturbative domain (Q = 2 GeV) at the expense 
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of parametrizing the effects of additional vector mesons by a simple factor K to be fitted 
to the data. 

We have presented self-contained formulae and parametrizations that allow the PDFs 
and F 2 of the photon to be evaluated as a function of a; and Q 2 , and also for non-vanishing 
target virtuality P 2 . A program with this information already encoded is obtainable on 
request from the authors. The anomalous PDFs are parameter-free; more precisely, the 
parametrization depends analytically on Q , A, and P 2 (besides, of course , Q 2 and x). 
Similarly, the state distributions can be evaluated for arbitray x, Q 2 , k 2 , and A. Based 
on fits to the available F 2 (x, Q 2 ) data, four different sets of hadronic PDFs are provided 
in order to illustrate the uncertainties in current PDF determinations. The sets differ in 
the value of Q (0.6 and 2 GeV) and the data included (in both cases only data above Ql 
are included in the fit), and the factorization scheme. The latter dependence is formally 
of NLO accuracy, but is found to be numerically significant in LO analyses. 

We have outlined the small-x behaviour of the PDFs and also discussed the description 
of the hadronic event properties in photon-induced reactions. A proper treatment requires 
knowledge beyond that of the usual inclusive PDF parametrizations. Moreover, we have 
investigated how the photon-to-qq splitting probability enters the eikonalization of jet 
cross sections. We have derived a formula that ensures that the various contributions 
are not counted twice. We have shown that it is mandatory to change from the inclusive 
PDFs to the state distributions. 

Finally we have investigated the limitations of the approximation of F 2 (x, Q 2 ) by the 
"FKP-plus-TPC/27" expression. We found that (with correct parameter choice) the FKP 
formula is a good parametrization of the anomalous part of F 2 (x, Q 2 ) at large x but fails 
for x < 0.5. Concerning the parametrization of the hadronic part, we demonstrated the 
effects of the necessary Q 2 evolution of the TPC/27 parametrization and, most impor- 
tantly, pointed out the strong correlation between the cutoff scale pro of the FKP formula 
and the (size and shape of the) hadronic part: if pro is varied only in the anomalous 
part, it introduces an artificial dependence on an unphysical parameter, which must not 
be there. 
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Figure 1: Comparison of the fits with the available |28], |29], | 
The full curve is for set SaS ID, dashed for SaS 1M, dash-dotted SaS 2D and dotted SaS 
2M. Data at nearby Q 2 values are shown in the same frame, and the theoretical curves 
are evaluated at some average Q 2 value of that range. The small vertical bars through the 
curves indicate the x value where W = 1 GeV. The behaviour to the right of this point 
should be taken as indicative only, with the possibility of large higher-twist contributions. 



Figure 2: Subdivision of the full F% parametrization by component, compared with data 
at Q 2 ~ 5.2 GeV 2 . The total F% is shown by the full curve. The lowest dashed curve gives 
the VMD contribution, and the next lowest the sum of VMD and anomalous ones. The 
third dashed curve, which coincides with the full curve for the DIS fits, gives the sum of 
VMD, anomalous and Beithe-Heitler terms. For the MS fits, the full curve additionally 
contains the contribution of the C 7 term. Note that this last term is negative at large x. 



Figure 3: Comparison of parametrizations of the u-quark and gluon PDFs at Q 2 = 4,40 
and 400 GeV 2 . Full curve is for set SaS ID, dashed for SaS 1M, dash-dotted SaS 2D 
and dotted SaS 2M, as in Fig. [l]. Note that a logarithmic y scale is used for the gluon 
distribution. For Q 2 = 400 GeV 2 also the x axis is logarithmic, so as better to show 
differences in the small- x behaviour. 



Figure 4: Target-mass (P 2 ) dependence of the u-quark and gluon PDFs, for Q 2 = 10 GeV 2 . 
From top to bottom: \fP 2 ~ = 0, 0.2 GeV, 0.7 GeV, and 1.4 GeV. The top two frames are 
for set SaS ID, and the bottom two for set SaS 2M. The remaining two look similar 
qualitatively. 



Figure 5: The hadronic part of F^x^Q 2 ) as obtained by evolving the TPC/27 
parametrization of K7. The full curve is the x dependence at the input scale Q 2 = 
0.71 GeV 2 , dashed curve at 4 GeV 2 , and dash-dotted curve at 100 GeV 2 . 



Figure 6: Comparison between the FKP and our parametrizations of the anomalous u- 
quark distribution as a function of x. The x and xu(x) values are shown in linear scale to 
the left and in logarithmic scale to the right. Common values are Q 2 = 4 GeV 2 , P 2 = 0, 
rif = 3, and A 3 = 230 MeV; additionally Qq = 0.36 GeV 2 in our parametrization. 

a) The full curve is our parametrization, dashed curve is FKP with i max = Q 2 /x and 
to = Ql/(l — x), and dotted curve is FKP with t max = Q 2 and to = Qq. 

b) The full curve is our parametrization plus our direct term Q58D, dashed curve is FKP 
with their direct term (|57|), m q /pTo = 0.3/0.52, t max = Q 2 /x and to = Qo/(l ~~ x )i an d 
dash-dotted curve is FKP as above except that to = Qo- 

c) The full curve is our parametrization plus our direct term (|58|), large-dotted curve is our 
parametrization plus the standard direct term (|59|), dotted curve is our parametrization 
plus the modified direct term (|60"D, dashed curve is FKP with their direct term fl5T|), 
m q/Pro = 0./0.6, t max = Q 2 /x and t = Ql/(l — x), and dash-dotted curve is FKP as 
above, except that t Q = Q 2 



o- 



26 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503384vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503384vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503384vl 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503384vl 



This figure "figl-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503384vl 



This figure "figl-6.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503384vl 



This figure "figl-7.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503384vl 



TPC, Q 2 = 0.7' GeV 2 



0.2 




0.2 



curves for Q 2 = 0.71 GeV 2 i 
i i i I i i i 0. 

0.5 1 



• PLGTO, 2 = 2.40 GeV 2 
□ TPC, Q 2 = 2.83 GeV 2 



JPG, d 2 = 1 .3T-.GeV 2 \ 




curves for Q 2 = 1 .31 GeV 2 \ 

I I I I I I I I L 

0.5 1 



0.25 




curves for Q 2 = 2.60 GeV 2 



0.5 



• TPC, Q 2 = 5.09 GeV 2 
□ TOPAZ, Q 2 = 5.1 GeV 2 
T PL0T0 , 2 = 5.30 GeV 2 



0.5 



0.5 



0.25 - 



• PLUTO, Q 2 = 4.30 GeV 2 




curves for Q 2 = 4.30 GeV 2 



0.5 




0.5 - 



• OPAL, Q 2 = 5.90 GeV 2 
□ CELLO, Q 2 = 7.00 GeV 2 




curves for Q = 6.50 GeV 



0.5 



Figure 1: Comparison of the fits with the available [28, 29, 30, 19, 31, 32, 33, 34] data. 
The full curve is for set SaS ID, dashed for SaS 1M, dash-dotted SaS 2D and dotted SaS 
2M. Data at nearby Q 2 values are shown in the same frame, and the theoretical curves 
are evaluated at some average Q 2 value of that range. The small vertical bars through the 
curves indicate the x value where W = 1 GeV. The behaviour to the right of this point 
should be taken as indicative only, with the possibility of large higher-twist contributions. 
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Figure 1: Continuation from previous page. 
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Figure 2: Subdivision of the full F^ parametrization by component, compared with data 
at Q 2 5.2 GeV 2 . The total F^ is shown by the full curve. The lowest dashed curve gives 
the VMD contribution, and the next lowest the sum of VMD and anomalous ones. The 
third dashed curve, which coincides with the full curve for the DIS fits, gives the sum of 
VMD, anomalous and Beithe-Heitler terms. For the MS fits, the full curve additionally 
contains the contribution of the C 1 term. Note that this last term is negative at large x. 
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Figure 3: Comparison of parametrizations of the u-quark and gluon PDFs at Q 2 = 4,40 
and 400 GeV 2 . Full curve is for set SaS ID, dashed for SaS 1M, dash-dotted SaS 2D 
and dotted SaS 2M, as in Fig. 1. Note that a logarithmic y scale is used for the gluon 
distribution. For Q 2 = 400 GeV 2 also the x axis is logarithmic, so as better to show 
differences in the small- a; behaviour. 
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Figure 4: Target-mass (P 2 ) dependence of the u-quark and gluon PDFs, for Q 2 = 10 GeV 2 . 
From top to bottom: \[P^ = 0, 0.2 GeV, 0.7 GeV, and 1.4 GeV. The top two frames are 
for set SaS ID, and the bottom two for set SaS 2M. The remaining two look similar 
qualitatively. 
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Figure 5: The hadronic part of F^x^Q 2 ) as obtained by evolving the TPC/27 
parametrization of i^. The full curve is the x dependence at the input scale Q 2 = 
0.71 GeV 2 , dashed curve at 4GeV 2 , and dash-dotted curve at 100 GeV 2 . 
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Figure 6: Comparison between the FKP and our parametrizations of the anomalous u- 
quark distribution as a function of x. The x and xu(x) values are shown in linear scale to 
the left and in logarithmic scale to the right. Common values are Q 2 = 4GeV 2 , P 2 = 0, 
rtj = 3, and A 3 = 230 MeV; additionally Ql = 0.36 GeV 2 in our parametrization. 

a) The full curve is our parametrization, dashed curve is FKP with £ max = Q 2 / x and 
= £?o/(l ~~ x )i an d dotted curve is FKP with £ max = Q 2 and t = Qo- 

b) The full curve is our parametrization plus our direct term (58), dashed curve is FKP 
with their direct term (57), m q /pr = 0.3/0.52, £ max = Q 2 / x and t = Qo/(^ ~ x )i an d 
dash-dotted curve is FKP as above except that £ = Qo- 

c) The full curve is our parametrization plus our direct term (58), large-dotted curve is our 
parametrization plus the standard direct term (59), dotted curve is our parametrization 
plus the modified direct term (60), dashed curve is FKP with their direct term (57), 
m<i/pTo = 0./0.6, £ max = Q 2 1 x and £ = Qo/(^ ~ x )i an d dash-dotted curve is FKP as 
above, except that £ = Qo- 



32 



